First-principles calculations and bias-dependent STM measurements at the 
a;-Sn/Ge(lll) surface: a clear indication for the 1U2D configuration 
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The nature of the Q:-Sn/Ge(lll) surface is still a matter of debate. In particular, two possible 
configurations have been proposed for the 3x3 ground state of this surface: one with two Sn adatoms 
in a lower position with respect to the third one (1U2D) and the other with opposite configuration 
(2U1D). By means of first-principles quasiparticle calculations we could simulate STM images as 
a function of bias voltage and compare them with STM experimental results at 78K, obtaining an 
unambiguous indication that the stable configuration for the Q:-Sn/Ge(lll) surface is the 1U2D. 
The possible inequivalence of the two down Sn adatoms is also discussed. 
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The Sn/Ge(lll) system has been thoroughly studied 
for its very rich phase diagram as a function of Sn thick- 
ness and sample thermal treatment [H, Q. In particular, 
the a-phase is obtained by evaporating 1/3 ML of Sn 
on a clean c(2x8) Ge(lll) surface after a short anneal- 
ing at about 500K. This surface is characterized by Sn 
adatoms regularly located on one out of three T4 sites 
of the bulk terminated Ge(lll) surface, resulting in a 
(-\/3x -\/3)^30° reconstruction. A particular interest for 
this phase arose after the discovery of a gradual and re- 
versible phase transition to a 3x3 reconstruction below 
220K 3]. This well known, but not yet completely un- 
derstood, surface transition was initially explained as a 
Charge Density Wave (CDW) formation below a critical 
temperature d, 0, H, 0]. However, due to the absence of 
nesting at the Fermi surface [1, 0] and to the presence of 
two Sn-4d core level components with a 2:1 relative inten- 
sity ratio both above and below the transition tempera- 
ture [^, 0|, this model has been questioned. Alternative 
models describing the -s/SxVS to 3x3 transition as an 
order-disorder transition have been put forward. Among 
them, the "dynamical fluctuation" model [l^, 11 1, sug- 
gests that the 3x3 reconstruction is the ground state and 
the \/3x \/3 reconstruction, observed at room tempera- 
ture, results from thermally activated rapid vertical oscil- 
lations of the Sn atoms. Below the critical temperature 
the adatoms fluctuation is frozen in the 3x3 periodicity 
in which Sn adatoms have different heights with respect 
to the underlying Ge substrate giving rise to a buckled 
surface. 

The exact structure of the 3x3 reconstruction is still a 
matter of debate. In fact, two possible conflgurations 
of this surface have been proposed, one with two Sn 



adatoms in a higher position with respect to the third one 
(two adatoms in Up position and one in Down position, 
2U1D for brevity hereafter) and the opposite configura- 
tion (1U2D). STM results alone cannot tell whether the 
surface configuration is 1U2D or 2U1D, because imag- 
ing empty or filled electronic states results in a hon- 
eycomb (an apparent 2U1D) image or a complemen- 
tary hexagonal (an apparent 1U2D) one, respectively 
d, [l2| • In order to understand which configuration de- 
scribes such a system, surface-sensitive structural tech- 



niques yjl, 1 14i. 11 a. 1 161.1 171 , Sn-4d core level photoemission 
spectroscopy [Sj, |9|, ll8| , non-contact AFM investigations 
[3], theoretical calculations lol 12 . 2ol . 21, 22, 231 were 
used, producing conflicting results. 

In this paper we present first-principles quasiparticle cal- 
culations that, combined with a bias-dependent STM 
study, clearly demonstrate a buckled surface structure for 
the Q:-Sn/Ge(lll) system at 78K with one out of three 
Sn adatom displaced upwards, solving the issue of the 
1U2D versus 2U1D configuration. The inequivalence of 
the two down Sn adatoms is also discussed. 
First-principles calculations have been carried out at a 
first stage at the DFT level (in a plane-wave pseudopo- 
tential implementation [i^l) and in the Local Density Ap- 
proximation (LDA). The Sn/Ge(lll) surface has been 
represented in a repeated slab geometry consisting of six 
Ge layers of 9 atoms each, saturated by H atoms on the 
bottom layer and with Sn adatoms on top [25| . Quasi- 
particle corrections to DFT-LDA eigenvalues have been 
afterwards calculated in the GW approximation [2^ to 
allow a closer comparison of theoretical and experimen- 
tal results. It has been shown 27| that LDA eigenval- 
ues give already a good description of surface states in- 
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troduced in the band gap by tin for Sn/Ge(lll). How- 
ever, in order to obtain an accurate picture of the whole 
metal/semiconductor interface, it is necessary to go be- 
yond DFT since electron screening is largely space de- 
pendent in such an inhomogeneous system and implies 
sizeable quasiparticle effects to be taken into account for 
those energy structures which are associated to germa- 
nium [28j. STM images have been simulated using the 
Tersoff-Hamann model [i^, as energy- integrated GW- 
corrected local density of states at a fixed height above 
the sample, using an average tip-sample distance of 5 A. 
Calculations were performed on the two possible 3x3 re- 
constructions starting from initial configurations having 
a vertical buckfing of 0.4 A for both 1U2D and 2U1D 
systems. In both cases, as expected [13, 12, [i^, the 
relaxation converged to the same energy minimum, cor- 
responding to a 1U2D model with a vertical buckling of 
0.36 A between up and down tin adatoms [13] • Hence, we 
also considered a metastable 2U1D model having a ver- 
tical buckling of 0.20 A (a value close to the structural 
results reported in Ref. [l7|) for comparing the resulting 
band structure and STM images with the ones derived 
from the stable 1U2D configuration. The GW-corrected 
surface band structures calculated along high-symmetry 
directions of the 3x3 SBZ are reported in Figure [T] for 
the 1U2D and the 2U1D configurations. The three sur- 
face bands are associated to the three Sn dangling bonds: 
the "up" adatoms are characterized by a filled dangling 
bond while the "down" adatoms have partially occupied 
ones (3ll . [3^ . This is seen, for example, from the pro- 
jected density of states (PDOS) at the three types of Sn 
adatoms reported in Figure [21 Consequently, the two 
upper bands in Figure [1^ are associated to 2D adatoms 
and the third band at lower energy to lU adatoms. Con- 
versely, the 2U1D band structure reported in Figure [Th 
shows that the two lower bands describe 2U adatoms 
and the third one located at higher energy describes D 
adatoms. Such PDOS curves suggest that STM images 
obtained with very small bias voltages should mainly 
show the two D adatoms in the 1U2D case and the two U 




in the 2U1D case, resulting in both cases in a honeycomb 
pattern. This hypothesis is confirmed by simulating STM 
images for the 1U2D and the 2U1D reconstructions (re- 
ported in Figure[3^ and[3j3, respectively) as a function of 
the bias voltage V. Indeed, in both cases the simulated 
images obtained at bias voltages lower than 0.2 V show 
a honeycomb pattern in both empty and filled states. 
Increasing the bias voltage, in the 1U2D case (Figure 
[3K), the filled states images gradually revert to the ex- 
pected hexagonal pattern, passing through an apparent 
\/3x-y/3 reconstruction at about 0.27 V, while the empty 
states images preserve the honeycomb pattern. As a re- 
sult, the calculated images between 0.3 V and 1.0 V show 
the well-known complementary honeycomb and hexago- 
nal patterns (for empty and filled states, respectively) 
reported in many papers. Interestingly enough, a further 
unexpected transition from honeycomb to hexagonal is 
observed in the empty states images at higher bias volt- 
age. As a result, both the filled and empty states simu- 
lated images at 2.0 and 2.3 V show a hexagonal pattern. 
The 2U1D simulated STM images reported in Figure [3)3 
show an opposite behavior: increasing the bias voltage 
above 0.2 V, the honeycomb to hexagonal transition oc- 
curs in the empty states images, crossing the apparently 
fiat reconstruction at about 0.28 V, while the honeycomb 
pattern is maintained in the filled states images. Further 
increasing the bias voltage, a new hexagonal to honey- 
comb transition is observed in the empty states series, 
resulting in a honeycomb pattern for both empty and 
filled states images. Summarizing the results obtained 
by analyzing the simulated STM images, we found that: 
a) at low bias voltage (i.e. lower than about 0.2 V) for 
both the 1U2D and the 2U1D models STM images with 
honeycomb pattern are predicted in both empty and filled 
states; b) at intermediate bias voltage the simulated STM 
images of the 1U2D system confirm previous results ob- 
tained at ± Q.55V I4I showing honeycomb and hexag- 
onal patterns for empty and filled states images, while 
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FIG. 1: (color online). GW-corrected surface band structure 
for the 1U2D (a) and 2U1D (b) configurations along high- 
symmetry directions of the 3x3 SBZ (shown in the inset). 



FIG. 2: (color online). Projected density of states at Up and 
Down tin adatoms for the 1U2D (a) and 2U1D (b) configura- 
tions. 




FIG. 3: (color online). Simulated 5x3 nm STM images for the 1U2D (panel a) and 2U1D (panel b) configurations and 
experimental 5x3 nm STM images (panel c, T = 78K, I = 20 pA) as a function of bias voltage. All the experimental images 
were obtained on the same area showing a Ge substitutional defect. ES= Empty States; FS — Filled States. 



the 2U1D results give the opposite condition; c) at high 
bias voltage, calculations predict that the STM images, 
for both empty and filled states, should provide a picture 
of the true surface reconstruction (i.e. hexagonal for the 
1U2D and honeycomb for the 2U1D). 
STM measurements were carried out at 78 K using a Low 
Temperature UHV-STM (Omicron LT-STM, base pres- 
sure 5T0~^^ mbar) and a tungsten tip cleaned in vac- 
uum by electron bombardment. Germanium substrates 
were cut from Ge(lll) n-type wafers. A clean c(2x8) 
Ge(lll) surface was obtained using a standard sputter- 
annealing procedure A nominal 1/3 ML Sn deposi- 
tion was performed at RT, followed by sample annealing 
at about 200°C. The formation of the a-phase \/3x\/3 
Sn/Ge(lll) was checked by LEED and STM. 
It is well known since 10 years by now 3] that in STM 
measurements empty (filled) states images display a hon- 
eycomb (hexagonal) pattern. However, such an undis- 
puted evidence has always been obtained, to our knowl- 
edge, using intermediate bias voltage (usually ± IV). 
Here, in order to verify the theoretical predictions and 
to discriminate between 1U2D and 2U1D configurations, 
we acquired a series of STM images on the same sample 
area ranging from ±0.1^ to ± 2.3T^, reported in Figure 
[3};. The evolution of such experimental STM images is 
strikingly similar to the one reported in the simulated 
STM images series for the 1U2D configuration (Figure 
[HK). The first predicted transition from honeycomb to 
hexagonal at low bias voltage in the filled states series 
is clearly visible in the first three images collected in the 



0.1 ~ 0.3 V range. In particular, the 0.2 V image shows 
an apparent VSx \/3 reconstruction, as predicted by the 
simulated STM image at 0.27 V. Furthermore, the sec- 
ond expected transition from honeycomb to hexagonal at 
high bias voltage in the empty states series was verified 
as well in the experimental empty states STM images. 
All the evidences described so far clearly indicate a re- 
markable match between simulations obtained from ab- 
initio calculations and experiments, pointing out a 1U2D 
configuration for the Q;-Sn/Ge(lll) surface. 
The only remaining puzzling apparent contradiction 
comes from the Sn-4d core level spectroscopy that shows 
an opposite ratio of the two components with respect to 
the expected one. A recent paper [l^ suggests an ex- 
planation for such a mismatch, proposing that the two 
down adatoms are not equivalent. In that case the ex- 
perimental core level spectra could be successfully fit- 
ted using three components in the deconvolution of the 
Sn-4d peaks. Actually, though the calculated height dif- 
ference of the two down adatoms is negligible (< 0.01 
A), the inequivalence of the two down Sn adatoms is 
confirmed by both ab-initio calculations and STM mea- 
surements. Indeed, the profiles (reported in Figure H]) of 
the 1U2D simulated images show that in some cases the 
two down adatoms have a different apparent height. In 
STM measurements the inequivalence becomes evident 
at high bias voltage in filled states images. However, our 
calculations show that such inequivalence is very small, 
as demonstrated by the theoretical energy difference be- 
tween the 4d core level peaks calculated for the two down 
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Sn adatoms (AE=0.01 eV, much smaller than reported 
in Ref. [3, AE=0.16 eV). 

In conclusion, by comparing first-principles quasiparti- 
cle calculations and experimental STM results at the a- 
Sn/Ge(lll) surface, we obtained a clear indication about 
the configuration of such a system. All the evidences, 
namely the fact that first-principles calculations predict 
a stable system for the 1U2D configuration only and the 
comparison between bias-dependent simulated and ex- 
perimental STM images, unambiguously indicate that 
the stable configuration for the a-Sn/Ge(lll) surface is 
the 1U2D. Furthermore, our data confirm the inequiva- 
lence of the two Sn adatoms suggested by Tejeda et al. 
[l8| . although with a much smaller difference. 
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FIG. 4: (color online). Profiles of the simulated STM images 
for the 1U2D configuration along the line shown in the inset. 
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